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ABSTRACT 

Star clusters are born in a highly compact configuration, typically with radii of 
less than about 1 pc roughly independently of mass. Since the star-formation 
efficiency is less than 50 per cent by observation and because the residual gas is 
removed from the embedded cluster, the cluster must expand. In the process of 
doing so it only retains a fraction fst of its stars. To date there are no observa - 



tional constrains for fst, although A^body calculations bv lKroupa et al 



( 20011) 



suggest it to be ab out 20-30 p e r cent for O r ion-ty pe clusters. Here we use the 



data compiled by 
Ae/Be s tars and by 



Testi et al 



study of 



(1997 



de Wit et al 



de Zeeuw et al 



1998, 



2004, 



19991 ) for clusters around young 



20051 ) around young O stars and the 



19991 ) of OB associations and combine these mea- 
surements with the expected number of stars in clusters with primary Ae/Be 
and O stars, respectively, using the empirical correlation betw een maximal- 



stellar-mass and star-cluster mass of 



Weidner fc Kroupa 



(|2006l ). We find that 



fst < 50 per cent with a decrease to higher cluster masses/more- massive pri- 
maries. The interpretation would be that cluster formation is very disruptive. 
It appears that clusters with a birth stellar mass in the range 10 to 10'^ Mq 
keep at most 50 per cent of their stars. 



Key words: stellar dynamics - stars: early-type - stars: formation - Galaxy: 
formation - open clusters and associations: general. 
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1 INTRODUCTION 

The birth sites of star clusters are dense molecular clouds. As the star formation efficiency, 
SFE = e = Meci/(Meci + Mgas), where M^d and Mgas is the mass in stars and gas, respec- 
tively, just before star formation ceases, is generally believed to be low, the expulsion of the 
surrounding gas leads to the release of a large fraction of the stars from the clusters. 
The most efficient mechanisms 



or gas expulsion are the radiation and winds of massive 
stars and supernovae. iKroupal (120051 ) gives a simple example which shows that the luminosity 
output of the OB stars should be strong enough to destroy their natal cloud, as e.g. a 15 
Mq star releases as much as 3 x 10^° erg per 0.1 Myr into its surrounding medium while a 
cluster of 10^ Mq has only a binding energy of 8.6 x 10^^ erg. Therefore, the cloud should be 
dispersed, even before the occurrence of the first supernova. An observational example is the 
Orion Nebula Cluster which is practically gas- free in its centre despite its young age of only 
about 1 Myr (iHillenbrand &: HartmannI 119981 ). But strong "luminosity leakages" through 



low-densit y holes in the rn olecular clouds may dampen the effect of radiation arising from 



OB stars (IDale et al 



20051). In this case, supernovae may play significant a role in dissolving 



molecular clouds (I Wheeler fc Bash 
Rapi d remova. 



1977 



Goodwin 



19971). 



of gas in t 



1978 



Hills 



le ca s e of e < 0.5 probably leads to th e total destruction of the 



1980 



clusters ( iTutukoy 
e ~ 0.2 — 0.4 (iNiirnberger et al. 



Mathieu 



2002 



1983 



Lada et al. 



Lada fc Lada 



3) 



1984). Observations show that 



20031 ) but bound star clusters like the 



Pleiades and the Hyades do exist. The problem of cluster survival after gas expulsion has 



been identified as one o 
With the use of the 



the key proble ms in astrophysical research ( jPavies et al.ll2006l ) 



AarsethI ( 1l999l ) A^bodyG algorithm and inclu ding 100 per ce r it pri - 



mordial binaries, stellar evolution and a realistic Galactic tidal field, iKroupa et al.l (120011 ) 
re-examined the problem of cluster survival after gas expulsion. They showed that from an 
embedded cluster with 10^ stars and brown dwarfs, a bound object with about 25 per cent 
of the initial number of stars can survive, even with e = 0.33 and explosive residual-gas re- 
moval. This remnant cluster is the core of an expanding OB association. Thus, the Pleiades 
and Hyades may have formed from a compact Orion-Nebula-type object but only retained 



St 



0.25 of their birth-stellar mass after removal of their residual gas. This fraction, 
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2nn3a 



b|), and also on the relative distribution of the stars and gas (lAdams 



cluster ( 


3oilv & KrouDa 


(Adams 


2000 


)• 



Kroupa fc Boilyl (12002) studied the effect of cluster dissolution on the observed initial 
cluster mass function (ICMF). They concluded that three different regimes of clusters may 
exist which they call type I, II and III clusters. 



(i) Type I. Sparse low-mass clusters which contain no O stars. They have stellar masses, 
Meci, below 300 Mq or N^d < 1000 stars. As the stellar winds and ionising radiation of 



the stars is low because of the lack of O stars, the gas ex pulsion time-sca. 



order as the crossing time of the cluster (a few Myr). Thus. IKroupa fc Boily 



e IS o: 



the same 



(12002h expected 



fst ~ 0.5 for type I clusters. 

(ii) Type II. Clusters with 10^ < N^^ < 10^ or 300 < Med ^ 30000 Mq. They have only a 
few O stars. But due to their still rather low mass, gas expulsion i s 'explosive' on a time - 



Kroupa fc Boilvl fl2002h 



scale of a few 10^ years. Given the destructive residual-gas expulsion 
suggested fst ~ 0.1 — 0.2 for type II clusters. 

(iii) Type III. Massive clusters with more than A^^eci ^ 10^, Med ^ 30000 M©. These clusters 
can have thousands of O stars but due to the high mass the ionised gas is expelled adiabat- 
ically, with a time-scale longer than the crossing-time of the cluster. If this is correct then 
fst ~ 0.5 for type III clusters. 



Kroupa fc Boilyl (|2002| ) faced was that fst was virtually unconstrained by 



The problem 
observations. 

The aim of this contribution is to show how the fraction of retained stars for type I 
and low-mass type II clusters may be observationally constrained in order to improve our 
understanding of their very early evolution and possible fate. To achieve this a large set of 



observations addressi ng the issue of c 



able in the literature (ITesti et al. 
19991) is used. 



1997 



ustering around young intermediate-rnass stars avail 



1998 



1999 



de Wit et al. 



2004 



2005 



de Zeeuw et al. 



In Section [2] we describe the procedure to calculate the surviving star fraction. This is 
followed by Section [3] where the results of this work are presented before they are discussed 
in Section m 
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2 THE PROCEDURE 



2.1 The stellar initial mass function 



The following multi-component power-law IMF is used to estimate the number of stars 
expected in a star cluster: 



-O!0 
-Oil 



= k{ 



( — )' 

i—Y 

(rno_\-"l ( my 
V mil / V "^0 / 

^ \ mil / V mQ J 

with exponents 



.miow ^ m < mil, 
, mil i^m < mo , 
, mo ^ m < mi , 
,7711 ^ ni < rrimax, 



(1) 



ao = +0.30 
Qi = +1.30 
02 = +2.35 
as = +2.35 



0.01 ^ m/MQ < 0.08, 
0.08 ^ m/MQ < 0.50, 
0.50 ^ m/MQ < 1.00, 
1.00 ^ m/MQ. 



(2) 



where dN = ^{m) dm is the number of stars in the mass interval m t o m + dm. The 



2001 



20021). The advantage 



exponents represent the standard or canonical IMF (iKroupal 
of such a multi-part power-law description are the easy integrability and, more importantly, 
that different parts of the IMF can be changed readily without affecting other parts. Note that 
this form is a two-part power-law in the stellar regime, and that brown dwarfs contribute 



about 4 per cent by mass only. A 



log-norrn a 



form below 1 Mq with a power-law extension 



to high masses was suggested by IChabrierl (120031 ). To day the observed IMF is und erstood 



to be an invariant Salpeter/Massey power-law slope (jSalpetei 



1955 



Massey 



2003D above 



; MfTi, being independent of the cluster densi 


ty anc 


metallicity 


or metaUicities 0.002 


(Massev & Hunte 


1998 


Sirianni et al. 


2000. 


2002; 


Pai 


■ker et al. 


2001: Massev 1998. 2002. 


2003 


Wvse et al. 


2002: ] 


Bell et al. 2003: 


Piskunov et al. 


2004: Pflamm-Altenbure & KrouDa 


2006 


• 



The basic assumption underlying our approach is the notion that all stars in e very cluster 
follow this same un iversal IMF, which is consistent with observational evidence (lElmegreen 



1999 



Kroupa 



200ll ). 



2.2 The maximum-star- mass vs cluster- mass relation 



In a series of rece nt publications ( jKroupa fc Weidner 



2003 



Weidner et al 



2004 



2004 



Weidner &: Kroupa 



20051 . 12OO6I ) the influence of stars forming predominately in star clusters which later dis- 



solve into the field was studied. During this process support for the possi bility of a maximum 
mass for stars was found on a statistical basis (for R136 in the LMC, IWeidner fc Kroupa 



Low-Mass Star Cluster Formation Efficiency 5 



as 
o 



0.5 
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_l I I L 



J I I l_ 
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Figure 1. The logarith m of the most-massive-star versus logarithm of the stellar cluster mass. The thick solid line is the 
semi-analytic result fromlWei dner &: Kroupal | |2004) . The thick dashed line is the result of the 'sorted sampling' Monte-Carlo 
experiment of|Weidner & Kroup? ('2006^. The dots with error bars are observational values for a sample of young clusters from 
the literature (sec Wcidncr & Kroupa 2006, for details on the list). The large trian gle is the most-massive star in a state-of- 
the-art hydrodynamical star cluster formation simulation llBonnell et al.ll2003l. liool) . The thin solid line on the left side of the 
corner, labelled with 'mmax = Med', indicates the limit where all mass of a star cluster is concentrated only in one star. 



2005 



2004 ) , a result la ter confirmed by several independent studies (iOey fc Clarke 



2005 



Koen 



Figer 



20061 ) ■ This work then yielded to a more thorough investigation of massive stars 



in star clusters, resulting in the finding of a probably physical (and not statistical) relation 
between the mass of the most massive star in a young ( < 3 Myr) star cluster and the mass 
of the harbouring star cluster (IWeidner fc Kroupal l2006l ) . One consequence of this relation 
would be an ordered formation of star clusters meaning that low-mass stars form first and 
that star-formation ceases with the appearance of the high-mass stars. This may be a natural 
outcome of termination of sta r forma t ion by feedback, and had been sugge sted in the study 



of the Hyades and Pleiades by lHerbigl (119621 ) and in a study of NGC 3293 by 
^19821 ). 



Herbst fc Miller 



The most-massive-star vs cluster-mass relation {thick solid line in Fig. [T]) then follows 
by using the cluster mass, Med, 



M 



eel 



m^{m) dm 



(3) 



miow 



and taking into account that there exists exactly one most massive star in each cluster. This 
condition can be written as 
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1 = / ^(m) dm. (4) 

Here miow = 0.01 Mq is the minimal fragmentation mass, mmax the mo st- massive star in a 
cluster and mmax* ^ 150 MfT^ the measured maxi mal stellar mass limit (IWeidner &: Kroupa 



2004 



Figer 



2005 



Oev fc Clarh 



2005 



Koen 



function 



m 



eel 



20061 ). On combining eqs. S] and [3] the analytical 

(5) 



is quantified by IWeidner Sz Kroupa! (12004J ) and shown as a thick solid line in Fig. [H 

With this relation we can now calculate the mass of the cluster for an observed most- 
massive (primary) star. 



2.3 The data 



Testi et al. 


(1997. 


1998. 


1999) 



main-sequence stars (Herbig Ae/Be stars) with a large set of near-infrared observations with 
two different methods. They looked for over-densities in the photometric i^'-band flux around 
these stars in comparison to the background away from the stars and directly counted the 
stars found in the i^-band close to these stars. Amongst their main results are: 

• Young star clusters are very compact (~ 0.2 pc). 

• The more massive a PMS star is the larger is the cluster around it. 

• There is a gradual change around the spectral type B7 from Ae stars without noticeable 
clusters to Be stars with clusters. 

Concerning the third result it is interesting to note that the transition mass is also roughly 



the sa me stellar mass above which standard accretion theory of stars starts to fail. iTesti et al 



(I1999I ) therefore conclude further that massive star formation "is influenced by dynamical 
interaction [s] in a yo ung cluster environ ment." 



The data set from 



Testi et al. 



(119991 ) already includes completeness limits down to which 
the masses of all stars should have been detected by their method. These limits are calculated 
for zero and 2 magnitudes of extinction in the i^'-band. Therefore the observed number of 
stars {Nohs in Tab.[l]) would be all stars above the first limit if no extinction would be present 
or Nohs would be only the stars above the second limit if there is extinction. But the data 
do not allow differentiation between these two cases. 



Testi et al. 



(119991 ) also already converted the colours of the A- and B-stars into effective 
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temperatures (Tefj) and lumi nosities. With the use o: 



mairi-sequence stellar models flD'Antona fc MazziteUi 



2001 



Meynet fc Maeder 



a large set of main-sequence and pre 



1997 



Hurley et al 



2000 



Behrend fc Maeder 



20031) the luminosi ties and Teg values are converted to stellar masses 



and ages here. Note that iTesti et al 
Mq and with the use of the 



1998r) h ave determ ined ages but only for stars below 6 



Palla fc Stahlej (11990 



19931 ) models. For a consistent approach 
all ages are re-derived here with newer models. Some of the new ages (~ 30% of the stars 
in the Testi sample) are substantiall y larger (rn o re tha n 200 %) than the older ones, but 



Testi et al. 



the bulk are reasonably close to the 
uncertainties in the ages a nd models fo r PMS stars. 

mm sample 



(1l998l ) values. This reflects the still large 



Complementary to the 



Testi et al 



de Wit et al 



(12004 I2OO5I) searched 



for evidence if O stars observed in the field originate from young star-forming regions or if 
in-situ field formation of massive stars is possible. Amongst their 43 candidate stars they 
found 5 with previously unknown small clusters surrounding them. These five stars are in- 
cluded in this study. The given sp e ctral types have been converted to Teg with the use of 



table 3 of 



Chlebowski &: Garmanyl (1l99ll ). The conversion of the lumin osities and 



into masses and ages are done wit h the same models a s applied on the 



Testi et al. 



sample. Of the remaining 38 stars 



de Wit et al. 



t he Tf,ff 
Jl998f ) 



(I2OO4J ) could trace back 27 stars to star 
forming regio ns - making them run-away stars which were dynamically ejected from young 



cluster cores (IRamspeck et al. 



2001 



remaining 11 stars is more puzzling. 



Pflamm-Altenburg fc Kroupal 120061 ). The nature of the 



de Wit et al 



(I2OOJ) classify them as O stars formed in 



isolation, representing the lower end of the cluster mass function. But they need a rather 
shallow cluster mass function with a slope of /5 = 1.7; a result in co ntrast to known ob 



servational values which are closer to ( 3 = 2 for various environments (ILada fc Lada 



Hunter et al. 



2003 



2003; 



Zhang fc Falllll999l ). Current star formation theories diffe r on the ques- 



tion i 



1998 



isolated forrnation of massive stars is possible - some argue again st it (IBonnell et al 



Stahler et al. 



agate the concept 



2OOOI: 



Bonnell fc Bate 



Yorke fc Sonnhalter 



2OO2I 



200 



3 



Ballv fc Ziiineckeiil2005l ) while others prop- 



Li et al 



2OO3I ). No unambiguous evidence for 



such star formation has been observationally found other than the existence of isolated O 
stars. 

Using OB associations with known open clusters as cores is even more difficult than the 
above described samples. As they are older (> 10 Myr) other dynamical effects (2-body 
relaxation, few-body encounters, tidal stripping due to the Galactic tidal field) already re- 
moved additional stars after all gas is lost from the remaining cluster. Furthermore, the most 
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massive stars may already have exploded as supernovae - making our method of estimating 
the cluster mass through the most-massive star impossible. Also, discrimination between 
members of the remaining open cluster and the OB association is difficult with current 
data. Additionally, an OB association can be the result of the dispersion of more than one 
cluster. This is actually seen e.g. in the Monoceros OBI association which harbours the open 
clusters NGC 2264 and Mon Rl and in Perseus 0B2 with the clusters NGC 1333 and IC 
348. Nonetheless, using the study of OB association s by 



de Zeeuw et al. 



Testi et al 



(119991 



19981) and the 



, an attempt 



de Wit et al 



i s mad e to compare to the results obtained from the 

(j2004j ) sample. Particularly suitable from that data set is the Cassiopeia- Taurus association 
with the open cluster a Persei at its core. For Cas-Tau 83 B stars are counted and for a 
Persei 30 B stars. This gives a total of 113 B stars of which 27% (/st = 0.27) are retained 
in the cl uster. An estimate o f the original cluster mass is more difficult. Using an age of 



50 Myr (Ide Zeeuw et al.lll999l ) for both cluster and association, the most-massive star still 



alive should be around 7.5 Mq according to the stellar evolution models described earlier in 
this work. The border between A and B stars is around 3.5 Mq. Assuming a canonical IMF 
(see Section [2?T|) . 0.68% of all the stars and brown dwarfs lie in the mass range between 3.5 
and 7.5 Mq. For 113 B stars this gives a total number of 16600 stars. With a mean stellar 
mass of 0.36 Mq for the canonical IMF, the resulting cluster mass is about 6000 Mq in 
stars. With the same method we analysed the assoc iation - cluster pairs Cepheus 0B2 - 



Trumpler 37 and Monoceros 0B2 - NGC 2244 from 



database and from 



de Zeeuw 



Massey et al. 



et al 



(1l999l ). supplemented 



(119951 ). All results are shown 



with data from the WEBDAH 
in Tabled! 

The relation between the maximum-star-mass and the cluster-mass from Section 12.21 
(eg. El) is then use d to connect the m a sses o f the most massive stars observed in the combined 



Testi et al. 



(Il999h and 



de Wit et al 



(120051 ) sample with the initial star-cluster masses. With 
the use of the canonical IMF (eqs. [1] and [2]) the number of stars expected for each cluster, 
Nexp, is derived for the two mass (completeness) limits, assuming 50% binaries in both cases. 
The stars in binaries are assumed to be chosen randomly from the IMF subject to the mass- 
constraint imposed by our relation in Fig. [H Table [1] shows the results of the conversion. 
Also shown in Table [J are the radii, rpH , of th e clusters as far as they have been deter- 



mined by 



Testi et al. 



mm or 



de Wit et al. 



(I2OO5I ) . In those cases where no radii were given 



^ http://www.univie.ac.at/webda// 
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(marked by " in Tab. [1]) the mean of the other radii has been chosen. These radii are needed 
to obtain the two-body relaxation times, trei, of the cluster with the use of the following 
formulae: 



^-1 = 0.1— [Myr] 



In(Ar) 



and 



2rcci 

0"ecl 



4 



/lOOMfT 



eel 



^ [Myr]. 



(6) 



(7) 



In Table [U we quote the estimated initial relaxation time, assuming the birth-cluster radii 
are r = 0.5 pc with Nf,xp,i stars, and the current relaxation time, assuming the measured 
cluster radii and the observed number of stars. As is evident, trei.in > age in the majority of 
the cases (66 %), such that the current low number of stars cannot be the result of evapora- 
tion from the cluster due to early two-body relaxation. However, trei,now < age indicates that 
the current remnant clusters are relaxation dominated. A^body modelling would be required 
to further quantify the sum of the effect of loss of stars through gas expulsion and through 
the later two-body relaxation after the remnant cluster has re-virialised. But since cluster 
disruption through two-body relaxation takes about 20 x trei,now in a solar-neighbourhood 
tidal field, the later relaxation-driven evaporation would n ot have had much time to act 



significantly. Stars without any evidence for clustering in the 



Testi et al 



(1l998l ) sample have 



not been included in Tab. [J as they are likely dynamically ejected from their birth place 



(jPflamm-Altenburg fc Kroupa 



Several of the stars of the 



2006). 



Testi et al. 



(119981 ) sample show considerable amounts of gas 
around them. All cases with an amount of gas larger than 50% of the initial embedded cluster 
mass (Meci) have been excluded from further analysis and are not listed in Tab. [U These 
clusters have not been evacuated from the gas yet and therefore have not lost stars through 
this process yet. Furthermore these large amounts of gas probably hide larger numbers of 
stars. In the remaining clusters only small quantities of gas are still left which may produce 
extinction but are not important for the dynamics of the cluster. 



3 RESULTS 



Kroupa fc Boilyi (120021 ) used a transformation function, /st,KB, to transform the mass func- 
tion of embedded clusters (ECMF) into the initial mass function of bound gas-free star 
clusters (ICMF). They used the following description: 
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Table 1. Numbe r, de signation, mas s limits with and without extinction, reel and the observed numbers of stars (Nahs) from 
iTesti et alj jl99^ and lde Wit et al.l ^QOS)- Ag e and star masses derived from stellar models, ted in ^.nd t^^i are calculated 
from eqs. |6]and[7] In the first case A^cxpl and r^^i = 0.5pc is used, in the second A^obs a^nd ^ccl- The initial embedded cluster 
masses (Med) are from the maximum-star-mass vs cluste r-mass relation (S ection [2. 2[l. while the expec ted number of stars 
(^exp li2) and the transformation factors (/^t 1,2) for the lTesti et alJ lll999l ) and the lde Wit et al.l ||2005|) sample are derived 
as described in Section 12.31 



Nr. 


Dcs. 


limit 1 


limit 2 


age 


mass 


Mod 


red 




trcl 


iVoba 


iVcxp.l 


iVcxp,2 




/at 






no ext. 


with ext. 








obs 


in 


now 








1 


2 






IMq] 


[Mq] 


[Myr] 


[Mq] 


IMq] 


[pc] 


[Myr] 


[Myr] 












1 


V645 Cyg 


< 0.1 ''^ 


0.5 


2.8 


29.2 


621.4 


0.6 


5.2 


0.2 


>5" 


588 


136 


0.0085 


0.037 


3 


MWC 137 


< O.l'' 


0.5 


3.3 


17.0 


268.5 


0.4 


4.1 


0.9 


>59'' 


265 


61 


0.22 


0.97 


5 


BHJ 71 


< o.l'' 


0.5 


0.4 


12.3 


165.4 


0.15 


3.6 


0.05 


4 


168 


39 


0.024 


0.10 


7 


AS 310 


< o.l'' 


0.5 


4.2 


17.5 


280.2 


0.4 


4.2 


0.6 


>37'' 


276 


64 


0.13 


0.58 


10 


BD-|-65° 1637 


< o.l'' 


0.5 


0.4 


7.5 


78.6 


0.4 


3.0 


1.0 


29 


84 


19 


0.35 


1.53 


11 


HD 216629 


< o.l'' 


0.5 


0.3 


9.2 


107.0 


0.1 


3.3 


0.1 


29 


122 


26 


0.26 


1.12 


13 


HD 37490 


< o.l'' 


0.5 


0.2 


10.1 


122.7 


0.14 


3.4 


0.1 


9 


127 


29 


0.071 


0.31 


18 


XY Per 


< 0.1 


< 0.1 


6.7 


2.9 


18.3 


0.08 


2.4 


0.06 


3 


23 


23 


0.13 


0.13 


19 


LkHa 25 


< 0.1 


0.16 


0.8 


4.6 


37.4 


0.3 


2.6 


0.5 


11 


43 


28 


0.26 


0.39 


22 


LkHci 257 


< 0.1 


0.29 


3.2 


2.9 


18.3 


0.3" 


2.4 


0.9 


15 


23 


9 


0.65 


1.67 


24 


VY Mon 


< 0.1 


< 0.1 


0.4 


5.4 


47.9 


0.25 


2.8 


0.6 


25 


54 


54 


0.46 


0.46 


26 


VV Scr 


< 0.1 


< 0.1 


2.2 


3.4 


23.4 


0.1 


2.5 


0.2 


24 


28 


28 


0.86 


0.86 
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6200 














0.25 




" For these clusters no 


Teei values 


were given 


in IXest 


i et alJ i 


ll99Sl) andlde Wit et 


al.l 1I2OO5I 


t). Therefore the 


mean r^^i 


of all 





clusters is used. 

F or these star s no m ass limits were given in lTesti et al. I (|l99l). 

In lTesti et al. 1 lll998h most of the mass limits were given as < 0.1 Mq. During the calculation of the expected number of stars 
these limits were ch osen as the h ydrog en burning limit of 0.08 Mq. 

For these clusters iTesti et al.l lll998l ) give only a minimum number of stars observed. Therefore the fst values are in these 



cases only lower limits. 

The clusters from the lde Wit et ah I ll2004l2005h sample are included twice as in the paper not numbers of stars but densities 
of stars with large error bars are given. First they are show n with the mean A^pbg and then with the maximum A^obs • 
^ The fst values for these clusters in OB associations from lde Zeeuw et al are determined in a different way than the 

rest. See text for details. 

/st,KB(Meel) = 0.5 - OAg{lM,,V, CTu/eeP ^M^^f ), 

with 

GiiM^ci; criM..,;iM::r') = [~^^) , (9) 

and the constants ctim^^i = 0.5 and IM^^^^ = logiQ^M^^^^) = 4.0. Thus, for example, at the 
mass 10"^ Mq, /st,KB = 0.1, ie a 10^ Mq cluster only retains 10% of its stellar population 



(8) 



according to this parametrisation. 
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Figure 2. Transformation factor as a function of em bedded cluster mass for the sources from the lTesti et al. I lll999h sample 
(dots with solid error bars), the Ide Wit et al.l ||2004| ') sample (boxes wit h dashed error bars) and the Ide^eeuw et al.l lll999l ') 
sample (triangles with long- dashed error bars). For the lTesti et al] l ll999l) and the Ide Wit et al.l (|2004|) sample the upper ends 
of the "error bars" are determined for th e cases with assumed extinction, while for the lower limit of the "error bars" no 
extinction was assumed. The errors for the Ide Zeeuw et al.l jigggl l sample are rather arbitrary - simply assuming a square-root 
(Poisson) error in the obser ved number of OB sta rs in the association-cluster pairs. The solid line shows the transformation 
function (eq. [SJ adopted bv lKroupa &: Boiivl ll2002h . 




With the values from Table [T] it is now possible to define a similar quantity: but not 
a transformation function but transformation factors, fst, for each observed s t ar clu ster. 



This is done by div iding the observed number of stars, Nohs, from 



Testi et ai- 



de Wit et al 



(119991) and 



(120041 ) by the calculated expected number of stars, N^^p- These values are listed 
in Table [H 

Fig. [2] shows the transformation func tion (eg. ED a s a so lid line and the derived transfor- 



mation factors for each clust er from the 
bars, from the Ide Wit et ai- 



de Zeeuw et al. 



Testi et al. 



(119991 ) sample as dots with solid error 



(120041 ) sample as boxes with dashed error bars and from the 
(Il999l ) sample as triangles with long-dashed error bars. For the observations 
the case without extinction is used as a lower limit of the error bars and the case with 
extinction as an upper limit. 

To reduce the scatter in the observations in Fig. [2] the lower and the upper limits on the 
fst values are c o mbine d separately into four mass bins. The res ults are sho w 'n in Fig. [3] for 



the 



Testi et al 



de Wit et al 



de Wit et al 



(119991) sample only, in Fig. H] for the combined 



2004 



Testi et a 



( 2004) sample, and in Fig. [5]for all three samples combined (iTesti et al. 



3i 



mm and the 



de Zeeuw et al. 



1999 : 



19991 ). The dots connected by a dashed line are the lower 



12 C. Weidner, P. Kroupa, D. Niirnberger and M. Sterzik 



0.8 



0.6 



0.4 





1 1 1 1 1 


1 1 1 1 1 


1 1 1 




< 


► -< 

. , 1 , , 


1 

N 

\ 

\ 

\ 

\ 

S 

S 

, , 1 . . 


> < 

, 1 , , 


t 



1 1.5 2 2.5 3 

log,o(MVMj 

Figure 3. Like Fig.[2]but with lTesti et al.l lll999l') data binned into four mass bins. The lower limit (no extinction corrections, 
dashed line connecting the dots) and the upper limit (with extinction corrections, solid line connecting the boxes) are treated 
separately. The error b ars arc the variance in the individual mass bins. The solid line shows the transformation function (eq. (Sjl 
adopted bv lKroupaX Boilv ( 2002, '). 



limits and the boxes connected by a solid line are the upper hmits. The error bars for the 
dots and boxes are the variances in each indivi dual bin. Here again th e solid line (without 
dots) is eq. [HI the transformation function from lKroupa fc Boilyl (120021 ). 

The differences between Fig. [3] and Fig. H] are rather negligible - indicating a good agree- 
ment between the two samples. Because there are substantially larger uncertainties in the 
de Zeeuw et al.l (119991 ) sample (ages, the mass of the most massive star at 50 Myr, the tran- 
sition mass between A and B stars, and the extrapolation from only 0.68 % of the stars to 
the total cluster population), it is not included in the further analysis. 

In Fig. m it can be seen that the lower limit of the binned fst factors lie somewhat below 
eq. [HI while the upper limits coincide with eq. [HI Therefore it might be plausible to reduce 
eq. [H]to 0.4 for Med < 1000 M©. But the large error bars in this investigation do not allow 
a further determination. 



In a very recent publication (IWang &: Looneyll2007l ) the amount of young stellar objects 



(YSO) around Herbig Ae/Be stars is further studied with the use of archival 2MASS and 
Spitzer IRAC data. While most of their targets are already included in this study and the 
rest is too obscured by dust, the following additional data for VY Mon and VV Ser have 
been extracted. While they observe 26 stars for VY Mon (Testi: 25) and 22 for VV Ser 
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Figure 4. Same as Fig. [3] but with the combined pTesti et al.l l|l999l ) and lde Wit et al.l ll2004h data binned. 
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Figure 5. Same as Fig. [3] but with the combined pTesti et al.l lll999h . lde"Wit et al.l ||2004) and Ide Zeeuw et al.l l ll999l) data 
binned. 

(Testi: 24), they also give number counts for YSOs in the whole observed field. In the case 
of VY Mon they find 42 YSOs which is quite close to the 54 stars expected to be released 
by the cluster through gas expulsion (see Tab. [T]). For VV Ser they find 148 YSOs. This 
is actually substantially higher than the expected number of 28, possibly due to other star 
forming activity in that region. 
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4 DISCUSSION AND CONCLUSIONS 



This contri b ution shows that w i th the use of available observational data (iTesti et al 



1998 



1999 



de Wit et al. 



2004 



1997 . 



20051 ) it is possible to constrain the star-loss due to gas 



expulsion in modest star clusters. The extracted upper and lower limits are shown in Fig. HI 
The transformation function, /st,KB; is consistent with the upper end of the data. Thus, 
clusters with initial masses in the range 10 - 10^ M© appear to retain about 50% of their stars, 
although the uncertainties are sufficiently large to allow even smaller retainment fractions 
fst- Therefore such clusters would appear, in the stages after gas expulsion, as expanding 
associations. 

In summary, cluster infant weight loss is generally rather h igh, all clusters w it h Mpr] ^10^ Mq 
losing ^ 50 per cent of their stars (Fig. H]). According to Kroupa fc Boily (j2002 ) infant 
weight loss may increase with cluster mass, but while being consistent with the present 
data this is not required by them, although the dashed constraints in Fig. H] may be seen as 
lending some support to this notion. Further observations of these and other star-forming 
regions are needed to address the fraction of retained stars with more confidence. 
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